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Abstract 

We study the nonlinear classical dynamics of an electron confined in a double dot potential and 
subjected to a spin-orbit coupling and a constant external magnetic field. It is shown that due to 
the spin orbit coupling, the energy can be transferred from the spin to the orbital motion. This 
naturally heats up the orbital motion which, due to the presence of the separatrix line in the 
phase space of the system, results in a motion of the electron between the dots. It is shown that 
depending on the strength of the spin orbit coupling and the energy of the system, the electronic 
orbital motion undergoes a transition from the regular to the chaotic regime. 
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I. INTRODUCTION 



The development of nanoscience, nanotechnology, and of the design of appropriate sys- 
tems for quantum information devices have triggered a large body of studies on the 
efficient and the controlled states preparation of systems at the nanoscale. In particular, 
magnetic nanostructures were demonstrated to be highly promising both from an applied {?]- 



10| and fundamental nonlinear physics ll|, |12J point of views. In this respect semiconductor 



quantum dots with a spin-orbit coupling offer are a good case study |13l-ll5|. A key element 
in these systems is the coupling of the electron momentum to its spin via the spin orbit (SO) 
coupling. A momentum dependent coupling offers a new way of manipulating the electron 
spin by influencing the electron momentum via an external periodic electric field. This is 
the idea of the electric-dipole spin resonance proposed by Rashba and Efros for the electrons 



confined in lateral systems on the spatial scale between 10 and 100 nm 



151 ] . However, a 



strong external electric field can considerab 



y influence the orbital dynamics and drive the 



system far beyond of the linear regime 16j. Usually nonlinearity implies a complicated 



behavior and difficulties may arise describing the dynamics 



16] . However, the non 



17 



inearity 



20| that 



may lead, on the other hand, to a variety of interestingand subtle phenomena 
have been discussed for a range of physical systems l|-|6| . Two aspects are of interest here: 
First, what defines the regimes where the systems exhibit strongly nonlinear behavior and 
possibly chaos, and secondly, whether these nonlinear it ies and the chaotic behavior can be 
utilized for application as lead to qualitatively new features in the dynamics. 

Here we address the first of these questions for the nonlinear electronic dynamics where 
the driving force is a constant external magnetic field. Thus, we consider an effect opposite to 
the electron-dipole spin resonance protocol offered by Rashba and Efros In particular, 
we will consider the dynamics of the electron with SO coupling term, confined in the quantum 
dot and subjected to the action external magnetic field. Via the external magnetic field, one 
can directly control the spin dynamics and therefore eventually influence the orbital motion. 
We will show that under different conditions, depending on the fields parameters and the 
strength of the SO coupling, different types of the dynamics can be realized. These various 
types of the dynamics will be linked to the structure of the phase space of the system. We 
will demonstrate that using driving external fields one can switch the chaotic behavior in 
the orbital motion on and off. 
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II. MODEL 



We consider a model system: An electron with mass m e is being confined to a double 
quantum well described by a potential of the form U(x) = Uq [—2 {x/d) 2 + (x/cf) 4 ]. Here 
Uq is the energy barrier, and 2d is the distance between the minima. We will consider the 
low temperature limit in which case the orbital dynamics is not sensitive to the thermally 
assisted tunneling. For the particular values of the parameters d ~ 100 nm and Uq ~ 20 
meV this imposes the following restriction on the temperature T < 100 K. For the typical 
for GaAs electron effective mass m e = 0.067m the semiclassical tunneling probability is 
small, namely exp — 8 ^ 2d ™e M° >~ io~ 4 . Therefore, a classical consideration is justified. To 
quantify the SO interaction we use the Dresselhaus type coupling term H so = apa x , where 
p is the electron momentum, and we assume H = 1. Consequently the Hamiltonian of the 
system, with the applied external magnetic B(t) being parallel to the z- axis, reads 



2 

H = ^- + U(x)+apa x + ^ B gB(t)^, (1) 



where \xb is the Bohr magneton, and g is the electron Lande factor. Introducing the charac- 



teristic maximum momentum p max = y/2m e Uo we can estimate the maximal precession rate 



due to the SO coupling term fi™ ax = 2a\ /, 2m e Uo. While a constant magnetic field of the am- 
plitude Bit) = B, induces a spin precession around z-axis with the frequency Qb — ^bqB. 



Therefore, if the amplitude of the magnetic field is strong enough B > 2ay/2m e Uo/ Hb9, 
Vt B > fi^ ax the spin dynamics is described by the relation a x (t) = cri°^ cos(f2^) and the 
Hamiltonian (OQ) takes on the form 
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H= — + U(x) +apcos(Q B t), (2) 
2m e 

(0) 



a x (t) = cos(n B *), = 1- 



In what follows, we will the treat time dependent term in (J2j) as a perturbation and for the 
sake of convenience we use dimensionless units by introducing the scaling: E — > E/AUq, 
x — > x/d, t — > Slst, p —> p/\ / 2mUo, a — > a^/m e /4f/ . 
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III. SOLUTION OF THE AUTONOMOUS SYSTEM: STRUCTURE OF THE 
PHASE SPACE 



Before treating the general time dependent problem, we consider the autonomous system 
and find the corresponding solutions. This allows to identify the topological structure of the 
phase trajectories. For the autonomous case a = the equation of motion corresponding to 
the Hamiltonian ([2]) reads 

P* = —fa = x ~ x ■ ( 3 ) 

Performing the integration in (j3J) and inverting the result, for energies E near the potential 
minimum we obtain an oscillatory behavior in the form 



x{t) = ±Jx\ + (x 2 - xf)cn 2 (x t/V2; kj, ~ < E < 0. (4) 

Here x ,i, x = y/l + y/l + 4E, X\ = \/l — y/l + AE, x\ = 2 — Xq are the turning points 
defined by the relation V(x) = E, cn(. . .) is the elliptic Jacobi cosine 2l|, and the parameter 
k has the form k 2 = (xq — x\) /x"q. Close to the separatrix we have E = and consequently 
k = 1. If the energy of the system is negative, then the electron is located in the left or in 
the right well and performs oscillation confined by the potential barrier. The oscillations in 
the different wells are described by the solutions with different signs in (j3J). The time of the 
oscillation is given by T = 2\^2K(k)/xo, where K{k) is the complete elliptic integral of the 



first kind 2l|]. If the energy of the system is positive, the solution of the equation ([3]) reads 



x{t) = x cn \^/x 2 - lt,l/kj , E > 0. (5) 

From the solutions (jl]), (J5]) we conclude that, depending on the values of the parameter K 
the dynamics has qualitatively different nature. They are separated by the value K — 1 
of the bifurcation parameter hinting on the presence of topologically distinct solutions (see 
FigJT]). In the case of a positive energy, the electron can overcome the potential barrier and 
undergoes inter-minima transitions, e.g. from the left to the right minima. The separatrix 
line that divides different types of solutions corresponds to the zero energy case E = and 
the separatrix solutions have the following form: 

cosh(t)' P cosh(t)' {} 
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FIG. 1. Phase portrait of the system ([3]). The phase trajectories with negative or positive values 
of the coordinate x < 0, x > describe the dynamics of the electron located on the left or the right 
well, respectively. These phase trajectories are separated by a separatrix line from the trajectories 
corresponding to the positive energy case. The separatrix crossing the nodal point (p = 0, x = 0) 
belongs to the area of a maximal chaos and arises if an external perturbation is applied. 

The dynamics near the separatrix is very sensitive to small perturbations. This fact can 
be exploited to help the electron overcomes the potential barrier and to perform a transition 
between two equilibrium minima. The role of such a perturbation in our case is played by 
the SO term in (j2J). The physical mechanism behind the switching of the electron position 
between the minima points is the formation of a homoclinic structure. Homoclinic structure 
is formed due to the applied perturbation in the vicinity of the separatrix crossing the nodal 
point. The formation of the homoclinic structure can be verified via the Melnikov function 



method 22J. The equations of motion for a perturbed system can be represented in the 



following general form 

• 9H °m f( A 

X = ~dp~ a f( x > p ' t )> 

9Hq . . , . 

Px = — -Q— + ag[x,p,t). (7) 

Here the terms H x = dHo/dx, H p = dHo/dp describe the unperturbed motion, while 
the contribution from the spin-orbit coupling are contained in the terms af(x,p,t) and 



ag(x,p,t). Taking into account ([2]) we immediately see / = cos(t) and g = 0. Using ([7]) we 
can write Melnikov's integral in the following form 

A(0)= J [H x (x,p)f(x,p,t) + H p (x,p)g(x,p,t)]dt. (8) 

— oo 

Melnikov integral is a measure of the distance between perturbed stable and perturbed 
unstable separatrix trajectories. Therefore, a change of the sign of the Melnikov integral is 
equivalent to the crossing of the stable and unstable separatrix trajectories. Consequently, 
the formation of the homoclinic structure can be identified as the change of the singe of 
the Melnikov integral. On the other hand, a crossing of the separatrix is equivalent to a 
transition between the minima of the potential well. Introducing a shift of the initial moment 
of time t — > t + 9 we can obtain the one dimensional parameterization for an ensemble of 
the separatrix trajectories and for the distances between the trajectories as well. Using the 
shifted separatrix solutions 

_ \/2 _ tanh(t + 0) 

X ~ cosh(t + 9) ' P ~ ~cosh{t + 9) , ( ' 

from (jBJ) we obtain: 

A(9) = , cos(fl). (10) 

v ; cosh(7r/2) y ' v ' 

From ( jTUl) we see that the dimensionless switching time is equal to the 9 ~ tt. Therefore, 
taking into account the connection between the real and the dimensionless time 9 — > 9 -Qg 
we conclude that the real switching time is proportional to the inverse precession rate 9 m 
tt/Qb ~ T. However, note that the switching happens only close to the separatrix and the 
system needs an additional time to reach the separatrix state. The time that the system 
needs to reach the separatrix state should be larger than 9 w tt/Qb- This can be checked 
by numerical calculations as well (see Fig. [2] and Fig. [3]) 

IV. WIDTH OF STOCHASTIC LAYER AND CRITERIA OF CHAOS 

As was mentioned above, the SO interaction should lead to the formation of a stochastic 



layer in the vicinity of the separatrix line 22j. The finite width of the stochastic layer 



can be evaluated using the theoretical approach developed in the papers 



23 



3- 



Since 



the separatrix divides the phase space in two parts a = (in, out), the perturbation V(t) 
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FIG. 2. Numerical solutions of the system described by eq. ([2]). The generated chaotic trajectory 
lies close to the separatrix E = 0.01, VLb = 1, oe — > ol/AUq = 0.2 
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FIG. 3. In order to reach the separatrix value of the energy E ~ 0, the system needs a time which 
is clearly larger than the inverse precession frequency t > 1/£Ib, fis = 1, a = 0.2 

apcos(flBt) leads to the formation of a stochastic layer on both sides of the separatrix, see 
Fig. HJ In what follows we refer to the inner and outward parts of the stochastic layer by the 
index (in, out), respectively. The inverse precession frequency 1/fls defines the time scale of 
the dimensionless time and therefore we set Qb — 1- 

The change of the energy of the system H Q and the angular variable if during one oscillation 
period T a forms the discrete map [Ei, tpi) —> (E i+ i m+1 ) , ti — > U + i = U + T a , a = (in, out). 

E i+1 = Ei + AH^( Vl ); Vi+1 = Vi + A^(E i+l ), (11) 
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Pcnurbed orbit 



FIG. 4. Perturbed trajectories in the vicinity of the separatrix [24j] 



where E { = H(x(ti),p(ti)), tpi = ip(x(ti),p(ti)) and AH*^, A(p a (E i+ i) are the increment of 
the energy and the phase during one period that may occur due to the applied perturbation. 
The change of the phase and the energy of the system Aip a (E i+ i) , AH a (ipi) can be evaluated 
via the following equations 

U+T a {Ei) 

A V ° ^Q B T a (E t+1 ), A£T(^)« J {H ,V}dt. (12) 



ti 



Here {H , V} is the Poisson bracket. Due to the nonlinearity, the frequency of the orbital 
motion depends on the energy of the system with w ff = 2ir/T a (E). Therefore, the resonance 
condition with the external perturbation in case of nonlinear system has the form m a bj a = 
n a flB- Here n a , m a are integer numbers. The fixed points for the particular resonance (m a : 
n a ) are defined via the following equations: 

E i+1 = Ei = E^ n % (13) 

A V a (E^ :n ^) « AH a ( V { ™° m ° ) ) = 0. (14) 

From (I13p one can define resonant values of the phase and energy E { r a:na) . Energy 

of the system itself can be presented as a sum of the resonant part ^ m<T, "" T ) an d deviation 
AEi 

Ei = El m ^ + AE^ (15) 



Using (IT51) we can rewrite map ( TTTj) in the form 

AE l+1 = AE t + AH^^), (16) 

dAcp(E i+1 ) 



(p i+ i = (fi + A(p a (E i+1 ) = ^ 



AE t 



„ _p(mo-:".CT) 
-C-i + l— -C<ct 



With the notations 

(17) 



_ _„(m CT :n CT )' 



Eq. (jTBl) can be rewritten in the following form 

I i+1 = Ii + Gjy^-AH'iipi), (18) 
Vi+i = fi + ^i+i- 

Using ([3]), (Q and (112]) one can derive the explicit expression for the change of the energy 
of the system Hq, during one oscillation period T a : 

U+T a {Ei) 

AH^tp,)^ [ {H ,V}dt= (19) 



k 



a 




2v / 2cosfi B t k f i 2\ / 2cos^l B t 

V2 / at — ^— — - 7-77 | = av 2 — — — — cos<^, 



cosh 3 [fi B (t - ti)} J cosh[n B (t - ti)] cosh(?r/2) 



In (fT9l) . due to the time localization of the profile of the solution (Q, we extended the 
limits of integration to infinity. For the evaluation of the phase increment Aip a (E i+ i) that 
occurs during one period, again we note that the separatrix divides the phase space on two 
parts cr = (in, out). For inner area a = in the period of the oscillation, for the solution 
(]!]) reads T = 2y/2K{k)/xQ. Therefore, for the phase increment during one period of the 
oscillation we have A(p(E a ) = To = 2\/2K{k)/xQ. Taking into account the logarithmic 
divergence of the elliptic integral K(k ~ In (16/(1 — k 2 )) /2 and the relations Xq = \/2, 
1 — k 2 = —Ei n = \Ei n \, the phase increment that occurs during one oscillation period, for 
the trajectory in the vicinity of the separatrix reads Aip(Ei n ) = Ti n = In (16/\Ei n \). Taking 
into account this fact, for the explicit form of the map (TTTj) we deduce 

^ +1 = ^ + ln(-^- ) , (20) 



E. 



E t+1 = Ei + aV2 — -rcosifi. (21) 

COSh(7T/2) 



In particular, for the fixed points of the first resonance (m a : l a ) from (1201) we deduce: 

E i+1 =Ei = E™ n \ <p i+1 = ¥ >i + 27im = ip™ 1 + 2nm. (22) 
From (1221) we obtain 

aV M^y = °' 1,1 im) = 2,m - (23) 

Taking into account ([2"3"]) we find 

<Pa = \ + to, k G N, \E%\ = 16e~ 2m \ (24) 

Using (p4|) we can rewrite the map (120|) in the vicinity of the fixed points in the following 
form 

= Ii- K m sin 6i, (25) 

Vi+l — + 

Here 

_ aV2 7T 

_ \Ej£\ cosh (7r/2) ' (26) 
is the coefficient of the stochasticity and for convenience we used the new variable with the 
shifted phase 6i — <fi + ir/2. The continuous limit of the map (123]) reads 

^L = -K m sm6, (27) 
di 

where the index i plays the role of time. The differential equations f )27|) can be derived 
straightforwardly from the effective Hamiltonian 

H cS = ^I 2 + K m sin6. (28) 

On the other hand, the effective Hamiltonian responsible for the map (26) reads 

H' eS = -I 2 - K mC osd W ~ v )- ( 29 ) 

v=— oo 

The separatrix energy for the Hamiltonian (12 81) is equal to the coefficient of the stochasticity 



for the map ([25]) E s = K m 22|. In addition, we note that the equivalence of our system 



(26) to the effective Hamiltonian (30) is very helpful. Since for the evolution of the width 
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FIG. 5. The phase portrait of the system given by eq. (I25p . About 200 trajectories are generated 
for the different values of the stochastic parameter K: a) K = 0.1, b) K = 0.7, c) K = 1, d) 
if = 1.2. 

of the stochastic layer for our system we simply can adopt the well-known results for the 
system 25J (30). In our derivations we have considered the trajectory that lies in the 
vicinity of the separatrix of the initial system (see FigJl]). Therefore, in our derivation 
E™ is supposed to be small however different from zero. The reason is that exactly on 
the separatrix E™ = of the initial system ([3]) a motion is impossible since the period of 
the oscillation diverges logarithmically T in = In (16/|-E^|) — > oo, E™ — > 0. Consequently, 
the coefficient of the stochasticity (|26j) is large but finite. Using (124)) we can rewrite the 
coefficient of the stochasticity in the following form 



K„ 



a 



V2 2) 



2mn 



7T 



16 



cosh (7r/2) 



(30) 
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From Eq. (!30|) we can easily infer the connection between the amplitude of the SO coupling 
a and the emergence of chaos. Chaos in the system appears if K m > 1. 

Taking into account (26) and (130]) we obtain 



a > 



4^o \ET n \ 
m e vrv/2 



cosh (tt/2) . 



(31) 



The width of the stochastic layer reads 25] 



5E m = 2[/ (47r) 4 v / i^exp 



7T 



(32) 



In ( 1HT|) and (13"2l we have used the definitions — >■ E/AU , a — > a^Jm e /4:UQ and from 
the dimensionless SO constant and the dimensionless energy we switched back to the real 
SO constant and the energy. The value of the width of the stochastic layer is a very im- 
portant quantity. Since SE m defines the width of the energy interval where chaos appears 
E > —5E m <. Note that the width of the stochastic layer 5E m depends on the two real 
physical parameters only: The barrier height Uq and the SO coupling constant a. Therefore, 
the energy when chaos appears in the system can be estimated easily and verified by the 
experiment. The width of the stochastic layer is proportional to the SO constant. Conse- 
quently, the energy window of the chaotic dynamics —5E m < E < is proportional to the SO 
constant as well. From (|3T|) we see that the chaos criteria K m > 1, connect to several real 
physical parameters such as the SO coupling constant a, the height of the potential barrier 
Uq, energy of the system, and the effective mass of the electron m e . Therefore, the effect 
of chaos should be easily observable in the experiment by tuning these parameters. Close 
to the separatrix when the energy of the system tends to zero E™ — > the chaos criteria 
K m > 1 holds even for the very weak SO interaction which confirms the sensitivity of the 
system near the separatrix. The simple relations (126]) -( 132]) naturally define a particular class 
of the materials where the transition from the regular to the irregular orbital dynamics of the 
electron can be observed easily on the experiment. A phase portrait for the orbital motion of 
the electron, for the different values of the coefficient of the stochasticity is shown in Fig]5j 
As we see from Fig. 5 with the increase of the coefficient K = K m , the system undergoes a 
transition from the regular motion Figj5^i) to the irregular one displayed in Figj5b). 
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V. CONCLUSIONS. 



In the present project we considered a simple model relevant for spintronics. The electron 
spin is coupled to the orbital motion via a SO coupling term and is confined in a double 
quantum dot. For the control of the spin dynamics, we considered a driving protocol based 
on an external magnetic field applied along z axis. We have shown that if the amplitude 
of the applied field is strong enough B > 2ay/2mUo/ fiB9, the spin dynamics is periodical 
in time. In particular, the spin rotates around z axis with a frequency — ^sgB and 
for the x projection of the spin we have a x (t) = cos(f2#t). Due to this the orbital 
dynamics can be reduced to an effective, time dependent, one dimensional Hamiltonian 
model (J2J). Different types of the dynamics can be realized. If the energy of the system is 
negative — 1/4 < E < 0, the electron is localized in the right or in the left well and performs 
oscillations bound by the potential barriers. However, via the spin orbit coupling channel, 
energy can be transferred from the spin system to the orbital motion. This naturally heat 
ups the orbital motion of the electron and due to the presence of a separatrix line in the 
phase space of the system ([3]), the electron may change from one potential well to the other. 
We call this effect a spin orbit coupling induced "tunneling" (even though this process is not 
a tunneling in the quantum mechanical sense) to stress the fact that the electron can perform 
inter- well transitions. We derived simple analytical equations (|26|) - (|32|) for the appearance 
of chaos and showed how this regime can be reached via tuning real physical parameters 
such as SO coupling constant a, the height of the potential barrier Uq , the energy of the 
system, and the effective mass of the electron m e . Therefore, the effect of chaos should be 
easily observable on the experiment. We also proved that close to the separatrix when the 
energy of the system tends to zero E™ — > the chaos criteria K m > 1 holds even for a very 
weak SO interaction, confirming by this the extreme sensitivity of the system near to the 
separatrix. In addition we derived expression for the width of the stochastic layer ( 1321) . We 
proved that the width of the separatrix is proportional to the SO constant and inversely 
proportional to the energy of the system. 
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